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Temperature of solid particles in colloidal plasmas in space 
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Abstract—A general formalism is presented for calculating the internal temperature of solid particles 
(“‘grains’’) in a colloidal plasma in space. Some illustrative calculations of the grain temperatures as a 
function of the plasma kinetic temperature, the plasma density and the degrees of dissociation and 
ionization are presented. 


INTRODUCTION 


Ir IS NOW BECOMING INCREASINGLY APPARENT that situations exist in space 
where solid particulate matter (grains) coexists with a plasma phase in a thermal 
steady state with the solid having an internal temperature JT, which is sig- 
nificantly lower than the kinetic temperature 7 (assumed the same for all plasma 
species) of the plasma phase. The prominent examples are the solid grains in the 
interstellar medium, in the interstellar clouds, and in many HII regions—the size 
range of such grains being ~ 0.01-1 wm (see Wickramasinghe and Nandy, 1972). 
In fact, a significant fraction of all matter in the observable universe may exist in 
this colloidal state. Some authors hold that this state of temperature dis- 
equilibrium is a physically reasonable state for much denser regions such as 
stellar atmospheres and envelopes, including the much-discussed protoplanetary 
solar envelope (Lehnert, 1970; Arrhenius and Alfvén, 1971; Arrhenius and De, 
1973; De, 1973; Alfvén and Arrhenius, 1976). Such regions have also been widely 
suggested to be the locales of condensation of cosmic solids (see reviews by 
Anders, 1971; Wickramasinghe and Nandy, 1972; see also Woolf, 1973; Herbig, 
1974). Thus a need seems to have arisen of developing the formalism defining the 
relationship between the temperatures T and 7, in the general space environ- 
ment. The most comprehensive work to date in this area is that of Lehnert 
(1970), who has investigated not only the question of the relationship between T 
and T,, but also some aspects of phase transition in colloidal plasmas in space. 
We present here what we consider to be a more refined version of Lehnert’s 
preliminary analysis of the relationship between T and Ty. 


ELECTRIC POTENTIAL OF GRAINS 


A grain in a plasma acquires a negative electric potential which, expressed in 
volts, has the same order of magnitude as the plasma temperature expressed in 
electron volts. This is because the electron flux to any surface element in a 
plasma initially far exceeds the ion flux (the electrons having much larger 
thermal velocity than the ions), so that there is an accumulation of electrons on 
the surface. This accumulation proceeds until the resulting potential on the 
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surface provides sufficient attraction to the ions and sufficient repulsion to the 
electrons so that the electron and the ion fluxes to the surface become equal and 
a steady state is attained. The value @ of this steady state potential is given by 
(Spitzer, 1969; De, 1974) 
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where e is the electronic charge, m, and m; are the electron and the ion masses, 
and k is the Boltzmann constant. For a hydrogen plasma, one finds = 2.5. 
However, photoionization of the grain by UV photons from the neighboring 
stars tends to reduce the negative potential of the grain and in some cases even 
to drive the potential positive. In general, this effect becomes important in a 
circumstellar medium only in plasmas with very low electron densities where the 
influx of plasma electrons to the grain surface becomes comparable to the efflux 
of the photoemitted electrons. Calculations of the grain potential including the 
effect of photoionization have been presented by Feuerbacher et al. (1973) who 
have also defined the conditions under which the potential becomes positive. 


THE AMBIENT RADIATION FIELD 


The radiation field in which a grain in a cosmic colloidal plasma is immersed 
derives principally from three sources: (i) the radiation from the proximal star or 
stars modified by the intervening medium and diluted by the inverse square law, 
(11) the radiation from the plasma itself—free-free, free-bound, and bound-bound 
radiations (see, e.g., Allen, 1973 for analytic expressions for these radiations), 
and (iui) the radiation from the surrounding grains. The net radiation field can be 
estimated as a combination of these three sources for a specific problem. For our 
general discussion, we shall characterize this field by an intensity function f,—a 
function of the wavelength A of radiation. 

In the following analysis we shall idealize the grains as spherical blackbodies, 
with a modification necessitated by the fact that for very small grains, the 
efficiency for emission and absorption of radiation is impaired for wavelengths 
greater than the dimension of the grains. A simple model for this efficiency is 
(e.g., Wickramasinghe, 1967) 


Q, (a) = 27a/d, 27Ta<A<« 
= 0<A <27a, (2) 


where a is the grain radius. The absorption law given by Eq. (2) seems to agree 
well with the observed wavelength-dependent extinction of starlight caused by 
the interstellar matter. A more realistic analysis should include the departures 
from the blackbody behavior arising from the particular grain material. 


© Lunar and Planetary Institute * Provided by the NASA Astrophysics Data System 


Temperature of solid particles in colloidal plasmas in space 81 


THERMAL STEADY STATE: GRAIN TEMPERATURES 


The temperature of a grain is obtained from the thermal steady’ state 
requirement that the net rate of heat loss of the grain equal its net rate of heat 
gain. 

The rate of thermal energy input to a grain from any plasma species is simply 
the impingement flux of that species multiplied by the average energy input to 
the grain from each impinging particle. The thermal steady state for the grains 
(assumed negatively charged) is then expressed by equating the net rates of heat 
loss to the net rate of heat gain of a grain as follows (De, 1973, 1974, 1976) 
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Here B, is the Planck function, c is the velocity of light, and N, and C, are the 
number density and thermal accommodation coefficient for the a component of 
the plasma, where a =m, a, i, and e refer to the molecules, atoms, ions, and 
electrons respectively. The relative abundances of these components in the 
plasma as a function of temperature are found by assuming dissociation and 
ionization equilibria. The temperature T; is the effective temperature with which 
the ions impinge on the grains. Because of the attractive force of the potential @ 
on the grains, this temperature is higher than T and is given by (De, 1976) 


LS TI teal (4) 





However, the effective impingement temperature of the electrons is the same as 
their actual temperature. 

The left hand side of Eq. (3) represents the heat loss of a grain (per unit 
surface area) by radiation. The first term on the right-hand side represents the 
heat absorbed from the ambient radiation field. The second through fifth terms 
represent the heat gained from impacts with molecules, atoms, ions, and elec- 
trons respectively. The last two terms represent the heat gained from association 
and recombination on the grain surface (y, and y; being the fractions of 
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incoming atoms and ions undergoing association and recombination on the grain; 
W, and W; are the dissociation and ionization energies; the factor 1/2 with y, 
arises from the fact that two atoms combine to form one molecule). The factors 
in curly brackets are the fluxes (per unit area of the grain surface) of the various 
components of the plasma. Note that for a plasma containing a number of 
chemical species, the second through last terms on the right-hand side should be 
summed over all such species. In certain situations it may be necessary to 
include in the analysis the heat of condensation, evaporation, or chemical 
reaction. This can be done easily by adding appropriately to the left- or the 
right-hand side of the equation the flux of the species involved multiplied by the 
relevant energy per particle of the species. However, in most cosmic situations 
(where condensable substances form only a trace of the plasma) any such 
contributions to the right-hand side of the equation will be insignificant com- 
pared to the heat input from hydrogen. 

In Eqs. (3) and (4) we have assumed that the grains are negatively charged. In 
cases where the grain potential @ is positive, the subscripts i and e in the fourth 
and fifth terms on the right-hand side of Eq. (3) should be interchanged, and the 
factor [1+] in the last term should be replaced by the factor exp(— yw). 
Equation (4) would then give the effective impingement temperature T, of the 
electrons. 

Let us next assume that the radiation field f, may be considered to be 
equivalent to a blackbody radiation field at a temperature To, 1.e., Ty 1s the steady 
state temperature that a grain would achieve in the same region of space if it 
were not being impacted on by the plasma particles. This temperature is thus 
defined by 


[- pacaaa =| BT IQ (ayaa. (5) 


Now, the maximum possible values of the temperatures Ty) and T, in any 
colloidal plasma would be about 2000°K—the highest permissible temperature 
for a condensed phase. Even for a radiation field characterized by such extreme 
temperature, relatively small amount of energy is distributed in the Planck 
spectrum shortward of A = a, when a is of the order of 1000 A or less—a range 
which roughly covers the majority of cosmic colloidal particles. In this range, 
little error is introduced by replacing the expression in Eq. (2) by Q,(a) = 27ra/A 
for all wavelengths from 0 to «. With this assumption we obtain a law 
analogous to the Stefan’s law (De, 1973) 


C ° mrs 
G |, BaTQ(a)da = o'TS, 

(6) 
G | Ba(TQ(a)da = o' Te, 


where o’ = 2.66 o a ergcm”’ sec '°K ~~, o being the Stefan—Boltzmann constant. 
The accommodation coefficients used in Eq. (3) have a somewhat wider 
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definitional scope than the conventional accommodation coefficients (e.g., 
Kaminsky, 1965) which are defined for small temperature differences between 
the gas and the solid phases. It appears, however, that accommodation 
coefficients are relatively insensitive to the gas temperature—at least from a 
theoretical point of view (McCarroll and Ehrlich, 1963). Thus in absence of any 
experimental information about the coefficients for the various plasma species at 
a given temperature, the best one could do is to use the theoretically calculated 
coefficients, or those determined experimentally for a certain temperature, for all 
plasma temperatures—until the temperature is so high that sputtering and related 
kinetic effects come into play. 

The main difference between the present analysis and Lehnert’s analysis 
(Lehnert, 1970) is that unlike Lehnert’s approximate expressions, we have 
calculated the heat flux from the plasma species to the grain surface exactly (De, 
1976). One of the resulting disagreements between the two analyses is that 
Lehnert has used a factor 3/2 in place of the factor 2 in the second through fifth 
terms on the right hand side of Eq. (3). The resulting error can in some cases be 
significant. 

It is easy to solve Eq. (3) numerically for the grain temperature using a 
simple iterative scheme. 


AN ILLUSTRATIVE CALCULATION 


We shall now make an illustrative calculation of the grain temperatures in the 
dense regions of the primordial circumsolar colloidal plasma envelope such as 
envisioned in the Alfvén—Arrhenius model (1976) of the origin of the solar 
system. Let us specifically consider the region between 2 and 4a.u. from the sun 
—which is of possible interest in the context of meteoritic science since this region 
is often suggested to be the locale of condensation of some meteoritic com- 
ponents. The colloidal plasma medium is considered to be transparent to the 
infrared wavelengths in which the grains emit and absorb most of their radiation. 
Thus, the radiation received by a given grain from the surrounding grains as well 
as from the plasma becomes a negligible heat source compared to the heat input 
from the impacting plasma species. Hence the temperature Ty owes itself to the 
radiation from the sun, and is taken here to be a typical 150°K. The grains are 
assumed here to have the size-independent limiting behavior of ideal black- 
bodies obeying Stefan’s law (i.e., Q, = 1 for all wavelengths). 

Since hydrogen is the dominant species in a cosmic plasma, we assume in our 
example that the plasma phase consists entirely of hydrogen. We assume the 
ionization of the gas to be due to electron collisions, and recombinations to be 
radiative. Under these conditions, ionization equilibrium is given by (Elwert, 
1952) 


——+— = 8.3 x 10 5 =P | - (7) 
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Dissociation equilibrium is given by the relation 


Na _ K(T) 


where the value of the dissociation constant K(T) as a function of temperature is 
given for instance by Tsuji (1964). 

The results of the calculation are shown in Fig. 1, where all accommodation 
coefficients and y, and y; have been set equal to 0.1. In this figure we have also 
plotted the degree of ionization as a function of the gas temperature. The gas 
density N labeled in the figure is the density if the gas were entirely molecular. 
We do not wish to emphasize the actual numbers in the figure, but the figure is 
illustrative in the sense that it demonstrates certain features of the plasma-grain 
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Fig. 1. Temperature of solid grains immersed in a tenuous hydrogen gas (plasma) in 
space as a function of the gas temperature and the gas density. The gas density N 
indicated in the diagram is the density if the gas were entirely molecular. With 
increasing gas temperature, progressive dissociation and ionization of the gas have 
been taken into account. The broken line represents the degree of ionization. The 
grains have been assumed to behave as ideal blackbodies. The ambient radiation field 
of the grains has been assumed to owe itself mainly to an external source (a proximal 
star) and to have such energy density that the grains would achieve a steady state 
temperature T,=150°K in absence of the gas. Other symbols in the diagram are 
explained in the text. 
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temperature relationship in the colloidal plasma considered. These features are: 


(1) For gas densities of the order of 10'°cm™ or less, the grain temperature is 
independent of the gas temperature in the entire range T = 2000- 
20,000°K. 

(2) The grain temperature is more sensitive to gas temperature in a partially 
ionized gas than in a neutral gas. 

(3) With N ~ 10% cm”, the grain temperature exceeds the melting point of 
even refractory substances for gas temperatures = 13,000°K. 

(4) The ‘“‘knees”’ of the curves in the gas temperature range 3000-6000°K 
arise due to the onset of dissociation resulting in substantial heating of 
the grains by the “atomic hydrogen torch” effect. A similar increase in 
heating occurs at the onset of ionization due to the large amounts of 
energy brought in by the ions or the electrons (depending on whether @ is 
negative or positive) and due to their recombination on the grain surface. 


REMARKS 


We have presented in this paper a generalized theory for calculating the 
temperature of solid particles in a cosmic colloidal plasma. We have shown that 
under certain conditions the solid grains can coexist with the plasma phase up to 
fairly high plasma densities and temperatures. The theory is also applicable to 
colloidal plasmas in other applications (e.g., Sodha and Guha, 1970) where, 
however, the plasma tends to be opaque so that the effect of the radiation field 
becomes significant. The result usually is that the difference between T and T, 
cannot be very large. 
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